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ABSTRACT 

We present a new analysis of the stellar population of the Orion Nebula Cluster (ONC) based 
on multi-band optical photometry and spectroscopy. We study the color-color diagrams in BVI, 
plus a narrow-band filter centered at 6200A, finding evidences that intrinsic color scales valid for 
main-sequence dwarfs are incompatible with the ONC in the M spectral type range, while a better 
agreement is found employing intrinsic colors derived from synthetic photometry, constraining the 
surface gravity value as predicted by a pre-main sequence isochrone. We refine these model colors 
even further, empirically, by comparison with a selected sample of ONC stars with no accretion and 
no extinction. We consider the stars with known spectral types from the literature, and extend 
this sample with the addition of 65 newly classified stars from slit spectroscopy and 182 M-type from 
narrow-band photometry; in this way we isolate a sample of about 1000 stars with known spectral type. 
We introduce a new method to self-consistently derive the stellar reddening and the optical excess 
due to accretion from the location of each star in the BVI color-color diagram. This enables us to 
accurately determine the extinction of the ONC members, together with an estimate of their accretion 
luminosities. We adopt a lower distance for the Orion Nebula than previously assumed, based on recent 
parallax measurements. With a careful choice also of the spectral type-temperature transformation, we 
produce the new Hertzsprung-Russell diagram of the ONC population, more populated than previous 
works. With respect to previous works, we find higher luminosity for late-type stars and a slightly 
lower luminosity for early types. We determine the age distribution of the population, peaking from 
2 to ~ 3 Myr depending on the model, a higher age than previously estimated. We study the 
distribution of the members in the mass-age plane, and find that taking into account selection effects 
due to incompleteness removes an apparent correlation between mass and age. We derive the IMF 
for low- and intermediate-mass members of the ONC, which turns out to be model-dependent, and 
shows a turn-over at M < 0.2 Mq. 

Subject headings: stars: formation — stars: pre-main sequence — stars: luminosity function, mass 
function — open clusters and associations: individual (Orion Nebula Cluster) 



1. INTRODUCTION 

The Orion Nebula (M42, NGC 1976), with its associ- 
ated young stellar cluster (Orion Nebula Cluster, ONC), 
is the nearest site of recent massive star formation and 
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is therefore regarded as a prototype of dense star form- 
in g regions. At a dist ance of ^ 400 pc from the Sun 
(Sandstrom et al.ll2007t iMenten et alll2007l: Iffirota et al.l 
l2007f ). the region is projected on the northwestern side 
of the Orion A giant molecular cloud, but contrary 
to the majority of s tar forming regions in the Galaxy 
(jLada fc Ladal 120031 ) . it is only slightly extinguished. 
This is due to the presence of the most massive members 
of the young stellar population, 6'^Ori-C in the Trapez- 
ium cluster in particular, whose ultra-violet emission has 
carved an ionized cavity in the molecular material which 
now engulfs most of the young cluster. The high extinc- 
tion provided by the background molecular cloud (up to 
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Ay — 100 mag, iBergin et alj (|1996[ )'). the position of 
this region at the anticenter of the Galaxy {I ~ 209°), 
and the high projected distance from the galactic plane 
(6 ~ —19°) conspire to provide a rather low foreground 
and background stellar contamination, enabling us to 
perform extinction-limited studies of a young, nearly co- 
eval, relatively rich and uncontaminated pre-main se- 
quence (PMS) cluster (1 - 3 x lO'^ stars aged ~ 1 - 3 Myr, 
IHillenbrand etani200l . 

Despite these advantages, observational studies of the 
ONC are not without difficulties. The nebulosity which 
surrounds the ONC population and the presence of cir- 
cumstellar material (disks and envelopes) cause highly 
non-uniform extinction among the members which must 
be corrected for each individual star. Furthermore, on- 
going accretion, variability due to star spots or stellar 
instabilities, influence the observed optical colors, which 
may vary with time. In principle, an observational effort 
including simultaneous multi-band photometry over a 
large wavelength range and adequate temporal sampling 
to account for variability effects may overcome these 
difficulties, disentangling the non-photospheric emission 
from complete spectral energy distributions (SEDs). 
Even in the presence of perfect measurements of pho- 
tometry, and with consideration of the above effects, un- 
certainties remain in the conversions between colors and 
magnitudes to temperatures and luminosities, which dif- 
fer from those of main sequence stars especially for late 
spectral types. 

The combination of these uncertainties make the con- 
struction of an accurate H-R diagram (HRD) based on 
photometry alone for the PMS population misleading. 
The construction of the HRD is more robust when pho- 
tometry is complemented with spectroscopy, since spec- 
tral types can be easily determined even though mod- 
erate extinction, and the accretion excess can be di- 
rectly identified from the line emissions or measurement 
at high dispersion of spectral veiling. Nonetheless, also 
in this case, the uncertainty in the spectral type-Tcff 
relation and, in addition, the difficulty of obtain spec- 
tra in a crowded field with bright, variable background 
emission limit the accuracy in deriving the stellar pa- 
rameters. The ONC has bee n repeatedly studied both 
at visible waveleng t hs fe.g.. iHerbig fc Terndrud IT986I: 
Prosser et ahl 119941: iHillenbrand l 119971: iRobberto et al.l 



2004') an d in th e ncar-IR (e.g 'Ali & Dctov ' '199E 
Hillcnbrand_et_aJ,l ri998; Hill enbrand fc. Carpenter 200C. 
Lucas & Roche 2000; Luhm an et al.ir2000l : iSlesnick etlll 
2004t iRobberto et al.. .2009, ) with increasing field cover- 



age and sensitivity, probing the cluster well into the sub- 
stellar regime. From the point of view of constructing 
the HR D of the clus t er, th e most complete work remains 
that of IHillenbrandl (|1997l ) (hereafter H97). Using com- 
piled plus new VI photometry and optical spectroscopy, 
H97 derived the extinction for ~ 1000 stellar (i.e., above 
the H-burning mass) members and used the dereddened 
photometry to place the stars in the HRD. This made 
it possible to derive the stellar ages and masses through 
the comparison with PM S evolutionary mode ls. 

In a previous paper (|Da Rio et all l2009bl hereafter 
Paper l) we have presented the optical photometry of 
the ONC obtained in UBVI, Ha and in the 6200l 
medium-band filters using the Wide Field Imager (WFI) 



at the 2.2m MPG/ESO telescope on La SiUa. In this 
paper we use these data to construct the HRD for the 
intermediate- (1 M0 ^ M < 8 M©) and low- mass 
(0.08 Mq < i\/ < 1 Mq) end of the stellar popula- 
tion. As in H97, we complement the optical photome- 
try with spectroscopy, extending the spectral catalog of 
that work with new spectral types from slit spectroscopy, 
and also from narrow-band photometry. This allows us 
to improve over H97 extending the completeness of the 
sample of members placed in the HRD. The use of multi- 
color photometry provides advantages in estimating the 
contributions of both accretion luminosity and redden- 
ing to the observed colors. However, this is possible only 
together with an independent estimate of T^s- For late- 
type stars it is not enough to rely on stellar photome- 
try, in particular on the V and / bands alone, due to 
the degeneracy between reddening and T^s (H97). In 
principle, multi-band photo metric observations can help 
in breaking the degeneracy (|Romamello et al.ll2002D . but 
the presence of non-photospheric emission, such as, for 
example, the U and B excesses due to accretion, limits 
the accuracy of the results. The most reliable method 
to break the temperature-reddening degeneracy remains 
to obtain an independent estimate of the effective tem- 
perature of the stars by means of spectroscopy. Thus 
we derive the stellar parameters only for the sub-sample 
with available T^s estimates. 

The paper is organized in the following way: in Section 
[2] we summarize the dataset used, which is a UBVI and 
TiO-band photometry from ground based observations; 
in Section [3] we study the color-color diagrams for all 
sources in the catalog in comparison with intrinsic colors 
computed from synthetic photometry; this allows us to 
test and correct empirically the models in order to derive 
the intrinsic, photospheric, colors valid for the ONC. In 
Section HI we isolate the sub-sample for which a measure- 
ment of Teff is available, either from optical spectroscopy 
or narrow-band photometry, and using this subsample, 
in Section [5] we derive the new HRD of the ONC, and 
discuss the age distribution, the mass function, and the 
correlation between ages and masses using two families 
of PMS evolutionary models. In section|6]we summarize 
our findings. 



2. THE PHOTOMETRY 

We utilize our multiband optical catalog presented in 
Paper I to construct the HRD of the ONC. The catalog 
consists of a nearly-simultaneous photometry over a field 
of view of ~ 34' x 34' which comprises most of the ONC. 
We detect 2612 point-like sources in the / band; 58%, 
43%, and 17% of them are also detected in V, B, and U, 
respectively. 1040 sources are identified in the Ha band, 
and in the TiO narrow- band at 6200A. 

The photometry is not color-corrected towards a stan- 
dard (e.g, Johnson-Cousins) system, in order to avoid an 
additional uncertainty introduced by these transforma- 
tions, which would be prominent and non-linear due to 
the large difference of the WFI filter bandpasses com- 
pared to the standard bands. Therefore our UBVI cat- 
alog follows the WFI instrumental photometric system, 
calibrated only for zero-point derived from comparison 
to standard field observations in order to conform to the 
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Table 1 

Effective wavelength, FWHM, E.W. and Vega 
zero-points for the WFI photometric bands used 



Band 




FWHM 


E.W. 


Vega Flux 




(A) 


(A) 


(A) 


erg s~^cm~^A~^ 


U 


3470 


668 


639 


3.40 • lO"'-* 


B 


4541 


1139 


1058 


5.88 ■ lO-'* 


V 


5375 


884 


791 


3.86 ■ lO-'* 


I 


8620 


1354 


1431 


9.50 ■ 10-10 


TiO 


6217 


193 


190 


2.47 ■ lO-'* 



VegaMag conventior0. 

In Table[T]we summarize the characteristics of the WFI 
bands we use. For each band the effective throughput is 
calculated multiplying the filter passband^ by the aver- 
age CCD quantum efficiencjQ (Q.E.). This is analogous 
to Table 1 in Paper I, which included, by way of illustra- 
tion, central wavelength and FWHM of the filters, but 
not accounting for the Q.E. The averaged flux of Vega 
for each band 

J fvcs.SxdX/ Js^dX (1) 

is presented in column 5 of Table [TJ where is the 
band throug hput and /vcga a calibrated spectrum of Vega 
(|BohliD|[2007i) . This quantity, multiplied by the equiva- 
lent width (column 3), provides the zero-points of the 
photometric systems. For each band, this is the flux in 
erg s"^ cm"^ from a zero magnitude source when the 
filter throughput is normalized to have 100% response at 
the peak. 

As described in Paper I, the 50% completeness limit 
of our photometric catalog is at ^ ~ 20.8. In the same 
work we have also presented a spectral classification of 
217 M-type stars, based on a photometric index derived 
from the V, I, and TiO magnitude. 

The main improvements upon the photometric cata- 
log of H97, which shares a similar field of view of our 
optical photometry, are a slightly deeper detection limit, 
a smaller pixel size, a higher effective wavelength of the 
WFI /-band than the Ic band, closer to the peak of the 
SED for the latest spectral types. Moreover, whereas the 
H97 photometry was obtained merging data from differ- 
ent instruments and epochs, our catalog is fully based 
on WFI observations carried out in a single night. In 
this way the risk of introducing an additional uncertainty 
from different systematic effects affecting parts of the 
stellar sample is removed. 

3. ANALYSIS OF THE PHOTOMETRIC COLORS OF THE 

ONC 

For a star of a given spectral type, the broad-band 
colors as well as the integrated fiux and bolometric cor- 
rection (BC), depend in general on the Toff, \ogg and 
[Fe/H]. For main sequence and post-main sequence stars 
the dependence has been investigated both throu gh em- 
pirical studies (jFloweii 119961: lAlonso et all I1999D , rely- 
ing on stellar templates, and through stellar atmosphere 

^ i.e. according to our photometric calibration, Vega has zero 
apparent magnitude in all the WFI bands. 

^ http:/ /www. eso.org/lasilla/instruments/wfi/inst/filters/ 
http:/ /www. eso.org/lasilla/instruments/wfi/inst/filters/WFI- 
Filters-2006-09-18.tar.gz 



modeling (e.g.. lAllard fc Hauschildtl (|1995[ )). For young 
PMS stars, unfortunately, an extensive empirical study of 
the intrinsic photospheric colors, and their possible vari- 
ation compared to MS dwarfs, has never been carried 
out. This is mostly due to the fact that the observed 
fluxes for this class of stars are strongly contaminated by 
excesses caused by accretion activity and circumstellar 
material; moreover young (few Myr old) PMS popula- 
tions are usually somewhat embedded into their parental 
cloud, with extinction differentially affecting the cluster 
members. As a consequence it is rather difficult to ac- 
curately measure the intrinsic photospheric emission of 
PMS stars. 

Further, the effects of non-uniform foreground extinc- 
tion, those from the amount and composition of any cir- 
cumstellar dust which affects both the total extinction 
and the reddening law, and those from the possible pres- 
ence of on-going mass accretion cause a combination of 
reddening and blueing effects beyond the uncertainties 
in the intrinsic colors. The dependence of observed col- 
ors on the various physical parameters involved is de- 
generate. For these reasons, optical photometry alone is 
insufficient to derive stellar parameters for young stars 
and additional constraints are needed, for instance an 
independent estimate of T^g for the star. In order to dis- 
entangle or at least constrain the various effects, a precise 
knowledge of the underlying photospheric emission - i.e., 
the intrinsic colors ~ of young stars is required; we focus 
our discussion here first. 

Our multi-band photometry on the ONC can inform 
us as to the behavior of the optical photospheric colors 
of young stars. The observed population spans a broad 
range in Tcffand also features a range of disk and accre- 
tion properties, including diskless and non-accreting ob- 
jects which can be used to infer intrinsic colors as a func- 
tion of Toff. In this section we compare intrinsic colors 
predicted using synthetic photometry (Section 13. ip . with 
the observed colors of stars in the ONC. This allows us 
to both test the accuracy of the model predictions, and 
to constrain empirically those valid for the ONC (Sec- 
tion |321)- Additionally, we use synthetic photometry to 
model the effects of extinction and mass accretion on 
the integrated colors. The results of this analysis will be 
used in Section 2] to derive the stellar parameters of the 
members, for the subsample of the catalog with available 
ToffCstimate. 

3.1. Synthetic photometry 

Given the wide range of stellar spectral types covered 
by our survey, spanning from early types to late-M stars, 
we must consider different families of synthetic spectra 
which cover different ranges of T^s. The selection of ad- 
equate model atmospheres is particularly critical for late 
type stars (M-type and later), whose SEDs are domi- 
nated by broad molecular absorption features and pos- 
sibly contain some dust in their atmospheres that can 
affect the observed spectra and colors. 

For intermediate tempe r ature s we used the NextGen 
models (Hauschildt et al.l Il999( ). This grid of syn- 
thetic spectra, computed in the range 3000 K<T'off< 
10000 K using the ste llar atmosphere code PHOENIX 
(|Hauschildt et al.lll997D is known to reproduce fairly well 
the optical colors of stars down to the K spectral class. 
For late-M stars these models are not adequate to re- 
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produce the observational data. iBaraffe et al.l (|1998[ ) 
show that they are unable to fit solar-metallicity V vs 
{V — I) color- magnitude diagrams for main sequence pop- 
ulations. The problems set in for {V — I) > 2, which 
corresponds to Tcff < 3600 K, and can be attributed to 
sho rtcomings in t he tre atment of molecules. In particu- 
lar, lAllard et all (|2000f ) illustrate that the mam reason 
for the failure of NextGen spectra in matching the opti- 
cal properties of late-type stars is the incompleteness of 
the opacities of TiO lines and water. In that work they 
present new families of synthetic spectra (AMES-MT 
models) with updated opacities; the optical colors com- 
puted from these mod els were t ested against the VRI 
photometric catalog of iLegget^ ()1992[ ). finding a fairly 
good match with the data down to the lowest tempera- 
tures. This suggests some confidence in the ability of the 
AMES models to reproduce the optical colors of main se- 
quence dwarfs, at least in the V — I range. We anticipate 
that in this section (see Figure [S]) we will further demon- 
strate the good agreement between the synthetic ZAMS 
from AMES models and empirical BVI color scales valid 
for MS dwarfs. Although it remains to be demonstrated 
that the AMES models are appropriate also for the lower 
gravities of the PMS stars, the match for the main se- 
quence gives us some confidence in the use of this grid of 
spectra. 

We therefore use the AMES-MT models for the low 
temperatures in our set of atmosphere models instead 
of the NextGen spectra; although this substitution is 
strictly required only for Tcff ^ 4000 K, we adopt it for the 
whole 2000 K< Tes < 5000 K range so to have a smooth 
connection with NextGen at intermediate tempe ratures. 
For T 'pff> 8000 K we use the Kurucz models ()Kuruczl 
119931) . Due to source saturation, we lack the most lumi- 
nous members of the ONC in our photometric catalog, 
and therefore our sample is not populated in the early 
type end of the population. The metallicity adopted is 
[M/H]= 0, the only value available for the AMES-MT 
grid, which is in good agreement with the observations 
of ONC members (e.g. D'Orazi et al. 2009). 

The values of log g span from 3.0 to 5.5 (in logarithmic 
scale with g in cgs units) with a step of 0.5 dex. During 
the PMS phase stellar contraction leads to a decrease 
of the stellar radius and therefore an increase in log g 
according to the relation: 



5.5 



log .9 = log 



M f R 



Mq \Rp 



L0 



•50 



(2) 



and 50 = 2.794 • 10'' g cm s ^. Examples o f these 
relations are shown in Figure [D for ISiess et al.l ([2000l ) 
isochrones from 0.5 Myr to the ZAMS. 

In general increasing ages correspond to increasing sur- 
face gravities and for the M-type stars (Teff<4000 K) 
the difference in \ogg can be more than two orders of 
magnitude during the PMS evolution, changing from 
log 5 ~ 3.5 for a star aged a few Myr to log 5 ~ 5.5 for 
a dwarf M-type stars. For early spectral types the dif- 
ference is lower, about 1 dex (log 5 ~ 3.5 to log 5 ~ 4.5). 
We therefore create, by interpolating within our grid of 
synthetic spectra, a family of spectra corresponding to 
each Siess isochrone, as a function of temperature. On 
this family we perform synthetic photometry to compute 
absolute magnitudes in the WFI instrumental photomet- 
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ric system, the standard of the photometric catalog pre- 
sented in Paper I. 

The absolute magnitude in a photometric band S\ of 
a star with a spectral energy distribution F\ , and stellar 
radius R is given by 



Af5, = -2.51og 



R 
lOpc 



dA-, 



ZP: 



Sx 



(3) 

where is a reference spectrum that gives a known 
apparent magnitude ZPs^ ; in the Vegamag photomet- 
ric system, which uses the flux of a Lyr as reference, 
/° = Fa.vega and the zero-points ^Ps^ vanish accord- 
ing to our definition of WFI Vegamag system (Paper 
I). We also assume an extinction Ax = 0, since we are 
interested in constructing a general observable grid of 
evolutionary models for unreddened objects. Eq. (j3|) can 
be rewritten then as: 



Ms, =-5 log 



^0 



43.2337- 5 log( ^ ) +B(Fa, 5a) (4) 



51og(^ 



S(Fa,5a) 



where 



B{Fx,Sx) = -2.5 log 



AFA^AdA 

> (5) 



j XFxycs&SxAX 



The latter term can be directly calculated having the syn- 
thetic spectrum for every Tcff point of a given isochrone, 
a calibrated Vega spectrum and the band profile of the 
WFI photometric filter. The Vega normalization is per- 
forme d using a recent reference spectrum of Vega ()Bohlinl 
12001 . 

In Figure [5] we present three observed color-color plots 
for the ONC in B, V, I and TiO bands, together with 
the intrinsic colors computed from synthetic photome- 
try using our grid of spectra selecting logg values from a 
ISiess et all ()2000D 2 Myr isochrone and for the zero age 
main-sequence (ZAMS) of the same family. The 2 Myr 
age is illustrative here, but consistent with the age of 
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Figure 2. Color-color diagrams for the ONC (circles). The solid 
lines are the loci of intrinsic photospheric emission computed by 
means of synthetic photometry, assuming a 2 Myr Siess isochrone 
(blue line) and the Siess ZAMS (red line) to define logg as a func- 
tion of T(,ft, for Ay = 0. The reddening vectors corresponding to 
Ay = 2 mag for two choices of Ry are reported as a reference. 



Table 2 

Transformation between reddening parameters. 



Rv H,(WFI) = ,j^4^ /j^(wFI) = ,j^^ 



2.0 


0.69 


1.69 


3.1 


1.06 


2.06 


4.3 


1.31 


2.31 


5.5 


1.49 


2.49 


6.0 


1.53 


2.53 



ONC derived from an HRD once Si ess et all (|2000D evo- 
lutionary models and a shorter distance for the ONC are 
used (though at least twice larger than th e value found 
by H97 using iD'Antona fc Mazzitelll' (1994) evolu tionary 
models) . As shown in Hillenbrand ct alj ([20081 ) , Siess 
models tend to predict higher ages, relative to other con- 
temporary model isochrones, for a given PMS population 
in the HRD. Also, using a shorter distance for the ONC 
increases the age significantly from that derived by H97, 
by - 50% 

In the same plot we show the extinction direction, for 
Ay = 2. We have assumed two choices of the extinc- 
tion law, the "typical" Galactic one, described by a red- 
dening parameter Ry = Ay /{Ab — Ay) = 3.1 and the 
"anom alous" one, considered more appropriate for the 
ONC (jJohnsonl 119671 : iCostero fc Peimbertill970l ). whose 
flatter variation with wavelength is described by a higher 
value Ry = 5.5. Since our photometry is expressed in 
the WFI instrumental photometric system, we computed 
the transformation of the reddening parameters R\ for 
the WFI bands as a function of Ry, which by defini- 
tion refers to the Johnson photometric standard. We 
performed the computation, for several choices of the 
standard reddening parameter Ry = Ay /E(B — V) (in 
Johns on photometric system) applying the lCardelli et ahl 
(|1989f) reddening curve on our reference synthetic spec- 
tra, and measuring the extinctions in both Johnson and 
WFI photometric systems. Examples of the transforma- 
tion between Ry (Johnson) and Ry , Ri (WFI) bands are 
listed in Table [2 

A broad range of optical colors is exhibited by the ONC 
stars in Figure [H Stars to the far right of the model 
ZAMS lines are heavily reddened stars, while stars to the 
far left, especially along B — V axes, can be explained 
in terms of highly accreting stars, as we demonstrate in 
Section [H Stars of intermediate color require detailed 
interpretation. It is important to recall from above that 
ZAMS colors are not necessarily applicable to the young 
stellar population of the ONC. According to the mod- 
els, the photospheric colors of pre-main sequence stars 
should vary significantly with age, in particular in BVI 
colors. From Figure[2]it is evident that the model ZAMS 
is systematically redder in B — V than much our data, 
in a way that is unexpected from a distribution of ac- 
cretion rates. The effect is significant for large colors 
{V - I> 2), corresponding to Toff < 4000 K, or M-type 
stars. Apparently the ONC color sequence is the locus 
seen just to the left of the model ZAMS. This suggests 
either a problem with the ZAMS relation or that the in- 
trinsic colors of young pre-main sequence stars are bluer 
than main-sequence dwarfs. We now explore whether 
the atmosphere models we have selected are appropriate 
to describe the intrinsic emission of first the ZAMS and 
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Figure 3. Comparison between the intrinsic colors derived by our 
synthetic photometry for a 2 Myr and the ZAMS (black lines) and 
the em piric al colors valid for MS d warfs from Kenvon & Hartmannl 
([1993) and IVandenBerg &: ClemI Il2003i '). in the Johnson-Cousins 
BVIc bands. The modeled ZAMS, compatible with the empirical 
counterpart, suggests a fairly good accuracy of the atmosphere 
models we used, at least for log g values of dwarfs. 

then the young ONC stars. 

We vahdate that the models can reproduce the ob- 
served colors of main sequence stars by showing in 
Figure [3] the 2Myr model and the ZAMS model com- 
puted in the Johnson- Cousins BVIc bands using our grid 
of synthetic spectra, as well as empirical colors valid 
of MS dwarfs, frorn . [Kenvon fc Hartmannl (1995) and 
IVandenBerg fc ClemI ([2003^. We test the correctness of 
the model ZAMS in the standard Johnson- Cousins bands 
since we don't have at disposal an empirical color scale 
for MS dwarfs in the WFI phtometric system. We stress 
that due to the differences between Johnson-Cousins and 
WFI systems, the predicted colors are slightly differ- 
ent than those presented in Figure [2^, in particular the 
{B — V) color is larger than that for our WFI bands; 
however the qualitative behavior remains unchanged: a 
young isochrone is significantly bluer than the ZAMS. 
On the other hand the model ZAMS is compatible with 
empirical color scales derived for the same photometric 
system, with offsets much smaller than the difference be- 
tween the 2 Myr locus and the modeled ZAMS. This 
confirms that the ZAMS computed by means if synthetic 
photometry on our grid of atmosphere model is fairly ac- 
curate. 

3.2. Constraining the intrinsic colors for ONC stars 

Returning now to the WFI photometric system and 
Figure [2l the question becomes what is the appropriate 
set of intrinsic photospheric colors to utilize. In gen- 
eral, the non- uniform extinction of the members displaces 
each star in the direction of a reddening vector, spreading 
the location of any hypothetical ONC isochrone in color. 
We assume that the reddening distribution is similar for 
members having different spectral types; this is a fairly 
good approximation but for the lowest intrinsic luminosi- 
ties for which high Ay are below the detection threshold. 
Thus, the fiducial points of the observed distributions 
for the ONC population are expected to be parallel to 
any appropriate pre-main sequence model computed for 
Ay = 0, and located at redder colors than the intrinsic 
photospheric locus (i.e. shifted by a non-negative amount 
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Figure 4. BV I color-color plot including only sources with small 
photometric error (< 0.1 mag), membership probability known and 
> 98%, no H«, excess, modest extinction from H97 (0.1 < Ay < 
0.5) (dots). The data have been de-reddened based on the ex- 
tinctions measured from H97. As in Figure (2] the red and blue 
lines represent respectively the synthetic ZAMS and 2 Myr model. 
The green line is our empirical calibration of the intrinsic colors, 
intermediate between the two (see text). 

of extinction). 

Although reddening from a 1 or 2 Myr isochrone could 
explain the observed displacement of the observed color 
locus redward of the model isochrones but blueward of 
the ZAMS, it is not adequate to simply adopt a model 
isochrone since one of the goals of this study is to assess 
the age and age spread for the ONC population. It is 
instructive to study the constraints available from our 
own photometry on the actual intrinsic BVI colors that 
are appropriate. We define an empirical intrinsic color 
isochrone in BVI for stars in the ONC by considering 
a set of stars that are free of both blue excess due to 
accretion and reddening due to extinction. 

From Figure |2^,c it is evident that the 2 Myr syn- 
thetic colors, when reddened, are more compatible with 
the observed sequence than the ZAMS locus. However, it 
is hard to provide a quantitative estimate of this level of 
agreement. This is because the observed sequence shows 
a large scattering, mainly due to the differential extinc- 
tion but, as we will demonstrate in Section |SJ also be- 
cause of broad-band excesses caused by stellar accretion. 
To simplify comparison of the model colors with our data, 
we isolate stars that: 

• are candidate non-accreting objects: we utilize the 
Ha excesses presented in Paper I and we select the 
sources with no excess. 

• have modest extinction: we consider the extinc- 
tions from H97, and select sources with < Ay < 
0.5 

• are known ONC members: we consider the mem- 
berships tabulated in H97, mainly collected from 
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IJones fc Walked ()1988D . and select stars with mem- 
bership probability m > 98% 

• have small photometric errors: we restrict errors to 
< 0.1 mag in all BVI bands. 

Restricted in this way, the sample satisfying all of the 
above criteria contains '^21 stars and includes the known 
ONC members for which the observed fluxes are closest 
to simple photospheres. By dereddening their colors sub- 
tracting the modest extinctions measured in II97, these 
stars should lie on our color-color diagram along a se- 
quence which represents the intrinsic colors valid for the 
ONC. This color locus is illustrated in FigurelH where we 
also overlay the predicted intrisic colors from synthetic 
photometry for the ZAMS and the 2 Myr isochrone. 
While the intrinsic colors of this sample of non-accreting 
ONC stars are bluer than those along the ZAMS, they 
are not as blue as the theoretically predicted colors for a 
1-2 Myr population. 

We thus confirm that the intrinsic colors for the ONC 
are not in agreement with the dwarf colors, but rather 
bluer as qualitatively predicted by the 2 Myr isochrone. 
However, small offsets are still present with respect to the 
young isochrone that are color-dependent. In particular, 
while the 2 Myr model seems appropriate for the coolest 
stars (~ 3000 K, corresponding to {V — T) ~ 4), the 
observed separation of the ONC color sequence from the 
ZAMS locus occurs at slightly larger colors {V — I >2) 
than predicted by the 2 Myr model. We chose therefore 
to use the empirical data to refine the model isochrone 
colors so as to best match the data of Figure S) To this 
purpose we consider the ZAMS intrinsic colors for Tcff > 
3700 K, the 2 Myr synthetic colors for Toff < 2900 K, and 
a linear interpolation between the two for 2900 K<Teff< 
3700 K. The result is shown also in Figure IH green line; 
it is evident how this empirical correction of the intrinsic 
colors now allows for a good match to the selected sample 
of non-accreting non-extinguished members of the ONC. 

We refer hereafter as the reference model these 
empirically-validated intrinsic BVI colors, which can be 
considered a fairly good approximation of the photo- 
spheric fluxes, as a function of Tcff, for the ONC. 

3.3. Bolometric corrections 

As described above, intrinsic colors for 2 Myr and 
ZAMS stars were derived by means of synthetic pho- 
tometry performed on a grid of spectra using the appro- 
priate log 5 vs. Toff. For the intermediate \ogg values 
between the two models that is needed for our empiri- 
cal isochrone, we adopt ZAMS values for Tcff> 3700 K, 
2 Myr for Toff< 2900 K, and a linear interpolation in 
between) to provide the empirically constrained set of 
model spectra that match the optical colors of the ONC. 

We can use these reference spectra to also compute the 
bolometric corrections (BC) valid for our filters without 
relying on relations taken from the literature that, be- 
sides requiring color transformations, may not fully ap- 
ply to pre-main-sequence stars. Clearly, based on our 
data, we cannot validate the accuracy of the BCs we de- 
rive from synthetic photometry. However, based on the 
evidences we have presented in Section and 13.11 and 13.21 
suggesting that our reference model is adequate to de- 
scribe the photospheric emission of the ONC stars at op- 



tical wavelengths, we are confident that the BCs derived 
on the same spectra are similarly appropriate. 

The bolometric corrections are defined as the difference 
between the bolometric magnitude of an object, found 
by integrating the SED across all wavelengths, and the 
magnitude in a given band: 



Jfx{T)-d\ I Jfxi(D)-dX 



(6) 

where C is a constant equal to the bolometric correc- 
tion for the Sun, BCq, in that particular photometric 
band. The values of C are irrelevant in our case, since 
for obtaining the total luminosities in terms of L/Lq (see 
Section [5] and Eq. [7]), BCq is subtracted and vanishes; 
thus we impose C = 0. The integrals in Eq. [6] have 
been computed for the V and / bands for the reference 
spectra, and for consistency we considered as a refer- 
ence solar spectrum a NextGen SED interpolated on the 
grid for solar values with Tcff=5780 K, \ogg = 4.43, and 
[M/H] = 0.0. 

The results are shown in Figure [5] as a function of 
temperature. For comparison, we also show these re- 
lations computed in the same way but using the surface- 
gravity values of the Siess ZAMS and the Siess 2 Myr 
isochrone. As for the colors, important differences arise 
for Toff< 4000 K, corresponding to M-type stars. In 
these regime, the difference between the three models is 
up to 0.15 mag for BC/ and up to 0.6 mag for BCv- For 
the derivation of the H-R diagram (see Section O we uti- 
lize only BC/, and the measured difference of 0.15 mag 
between the 2 Myr model and the ZAMS translates into 
an offset of 0.06 dex in logL, amount small enough to 
be irrelevant for the computation of the stellar parame- 
ters. Therefore, we assume the reference model BCs for 
consistency with the model we use to derive the intrinsic 
colors, confident that possible uncertainties in the BCs 
do not affect significantly the result. The major advan- 
tage of our approach is that we can compute consistently 
these quantities for the WFI photometric system which 
our photometry follows. 

3.4. Extinction law 

We are making an implicit assumption that the extinc- 
tion in the entire region under study follows the same 
reddening law. Differences in grain size and composi- 
tion might be expected, due to the complex structure 
of the ISM within the Orion Nebula, the non-uniformity 
of the foreground lid of neutral material and, for a frac- 
tion of members, the presence of circumstellar disks in 
different evolutionary state. Optical photometry alone 
cannot provide an independent assessment of the most 
adequate reddening law for each source, given the degen- 
eracy for late-type stars between the observed colors and 
Teff and the presence of non-photospheric effects. From 
our color-color plots, however, we can make some general 
considerations. Whereas the {V — I) vs. {TiO — I) plot 
of Figure [5] does not provide useful information, as the 
reddening vectors for the two choices of Ry are nearly 
parallel due to the limited wavelength range spanned, in 
the {V — T) vs. {B — I) plot the two reddening vectors 
are clearly separated. For spectral types K and earlier 
{{V — I) < 2), our template model is nearly parallel to 
the two reddening vectors (i.e., the above mentioned de- 
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Figure 5. Upper panel: intrinsic color (V — I)o as a function of 
temperature, computed for our 2 Myr isochrone, the ZAMS, and 
our empirically calibrated reference model, in the WFI instrumen- 
tal bands. Middle panel and lower panels: Bolometric corrections 
derived in the same way, for the I and V WFI instrumental bands. 
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Figure 6. (B - V) vs (V - I) plot of the ONC population (dots), 
with the 2 Myr reference model (solid line). The two thick lines 
correspond to the reddening direction assuming Ry = 5.5 (upper 
line) and Ry = 3.1 (lower line) extinction laws, applied to the 
tangent points of the model. 



generacy between spectral types and reddening), mean- 
ing that all the highly extinguished members within this 
range should be located in a narrow stripe at higher color 
terms. 

In Figure [S] we show the reddening directions for the 
two values of Ry applied to the tangent points with the 
reference model. If Ry ~ 5.5 were the correct choice, 
no stars would be detected to the right from the associ- 
ated reddening direction, in contrast with our data. On 
the contrary, the Ry =3.1 reddening limit encompasses 
almost the entire photometry, with the exception of few 
scattered stars with high photometric errors in (B — V). 
Clearly, the scatter in our data and the photometric er- 



rors limit a precise estimate of the actual value of Ry. 
Nevertheless, we can conclude in an unquestionable way 
that the typic al galactic law is more compatible with our 
data than the iCostero fc PeimbertI ([l97()) . 

Given the uncertainties, and for the purpose of dis- 
cussing how the choice of Ry affects the scientific results, 
we perform the analysis presented in the next sections 
using both values, but we present mostly the results for 
Ry = 3.1, mentioning any systematic changes induced 
by assuming the higher where relevant. 

3.5. Mass accretion 

The color-color diagrams of Figure [5] show that there 
is a fraction of stars lying significantly to the left of the 
models. In this region, according to our reference model, 
a star of given temperature and positive reddening can- 
not end up. The fact that the TiO-I vs V-I frame is the 
"tightest" , while the [TiO-I) and {V-I) vs. [B-V) are 
more scattered suggests that these sources have signifi- 
cant excess in the B band. Moreover, low {V — I) stars 
are scattered only bluewards with respect to the refer- 
ence model, confirming the presence of a i?-band excess 
for a subsample of sources with low reddening. 

We investigated the possibility that the blue excess 
is due to accretion luminosity. According to the cur- 
rent paradigm for classical T Tauri stars (CTTS), gas 
is accreted onto the central star from an accretion disk 
truncated at a few stellar radii by the stellar mag- 
netic field. From the inner edge of the disk, the ma- 
terial fall s along infall columns l3eanied by magnetic 
field li nes (IKoenigllll99il:IShu et aLlfTOOllMuzerolle et al.l 
Il998bf) . Several lines of observational evidence sup- 
port the magnetospheric model, including: (1) the 
broad emission lines emitted by the infalling gas, pre- 
senting occasionally redsh i fted absorption components 
(|Calvet fc HartmannI IT99I iMuzerolle et a\\ \\^9Sa): (2) 
the near-infrared spectral energy distributions (SEDs), 
consistent with thos e of disks truncated at small radii 
(|Kenvon et al.l 119941: iMever et al] I1997D : (3) the vari- 
ations in luminosity consistent with the presence of 
hot spots, originatin g where the accret ion flow i mpacts 
the s tellar surface (jHerbst et al.l 119941 : iGuUbring et al.l 
|2000[) . The presence of ongoing accretion processes in 
th e ONC has b e en inv estig ated by several s tudies , e.g., 
in iRebull et all (|2000l ) and iRobberto et"all ([200l from 
the [/-band excess i s used t o obtain mass accretion rates, 
or in i Sicilia-Aguila r et al.l ([2(305) and Hillenbrand et al] 
(|1998| ) using optical spectroscopy. 

In order to verify the possible role of accretion also in 
our BVI bands, we simulated a typical accretion spec- 
trum to b e added to our re ference models. According to 
iCalv et fc'GuUbring (Il998| ). the total accretion column 
emission can be well approximated by the superposition 
of an optically thick emission from the heated photo- 
sphere below the shock, and an optically thin emission 
generated in the infalling flow. The SED of the opti- 
cally thick emission is compatible with a photosphere 
with Teff = 6000-8000 K, while the optically thin prc- 
shock column can be modeled by thermal recombination 
emission. The relative fraction of the two components 
with respect to the total emission produced by accre- 
tion is about 3/4 and 1/4, respectively. We therefore as- 
sumed for the optically thick emission a black body with 
T = 7000 K, whereas for the thin emission we used the 
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Figure 7. Accretion spectrum simulated with Cloudy. The sohd 
line is the total emission, which is the superposition of an opti- 
cally thick emission, with refi=7000 K, of the heated photosphere 
(dashed line) and the optically thin emission of ionized gas with 
density n = 10* cm~^ (dotted line). 

ve rsion 07.02.01 of the Cloudy software, last described 
bv lFerland et al.l ( )1998| ). to reproduce the spectrum of an 
optically thin slab (n = 10^ cm'^) also at T ~ 7,000 K. 
The derived, low resolution, accretion spectrum is shown 
in Figure [T] 

We calculated the displacement in the color-color di- 
agrams caused by the presence of ongoing accretion 
by adding these spectra, with an increasing ratio of 
iaccr/^tot, to our reference model. For any given com- 
bination, we recalculated the synthetic photometry, de- 
riving the {B — V), and {V — I) colors in the WFI in- 
strumental photometric system. The results (Figure [5]) 
show the derived loci for each stellar temperature, all 
converging to the point representative of the pure accre- 
tion spectrum, as well as for fixed iaccr/^tot ratios. 

The blue excess that is observed for a number of the 
ONC targets in the {V — I) vs. {B — V) diagram appar- 
ently can be explained by an accretion spectrum added 
to the stellar emission, with typical luminosity well below 
10% of the total luminosity of the star. Figure [5] clearly 
shows that small fractions of accretion luminosity have a 
larger effect on the colors of low-temperature stars, indi- 
cating that accretion effects must be taken into account 
especially in the study of the low-mass end of the ONC. 

We have tested how critical is the assumption of a 
unique accretion spectrum for the modeling of the color 
excess. When the accretion temperature is set to be 
9000 K instead of 7000K (see Figure El, the colors of the 
pure accretion spectrum (i.e. the location, in the diagram 
of Figure |5] where all the curves join) decrease by 0.15 
mag in {B — V) and 0.2 mag in {V — I). This is a small 
difference already in the limiting case of -/jaccr/-^tot=l, 
for much lower accretion, which is the case of the major- 
ity of our stars, the position of the accretion curves of 
Figure [5] remains almost unchanged. We conclude that 
our approximation of applying a unique accretion spec- 
trum to model the stellar colors contamination is enough 
robust for our purposes. 

4. ANALYSIS OF STARS WITH KNOWN SPECTRAL TYPE 

In the previous sections we have isolated the intrin- 
sic color scale which best fits the photospheric emission 
of the ONC members. From the intrinsic colors, taking 
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Figure 8. BVI color-color diagram of our WFI catalog. The 
thick line represents the synthe tic C MP obtained from the ref- 
erence spectra selected (see sub il3.2| l for Ay = and the perti- 
nent effective temperatures are shown along the curve. Thin blue 
lines represent the computed displacement, for a star of a given 
temperature, adding a component of the accretion spectrum (as 
in Figure O and increasing the parameter Laccr/itot from to 
1. Loci for Laccr/itot = 0.5, 0.1, 0.01 and 0.001 are overplotted 
(dashed lines). Reddening vectors corresponding to an extinction 
of Av = 2 are shown for the two choices of reddening parameter 
Ry assumed in this work. 

also into account possible color anomalies through the 
addition of an accretion spectrum, we can determine the 
extinction Ay and therefore the stellar luminosity. The 
combination of spectroscopic information with accurate 
photometry, together with our determination of the in- 
trinsic relations (intrinsic colors and BCs for the young 
ONC, and the treatment of Laccr excess), allow us to 
build the H-R diagram for the ONC improving over the 
analysis of previous works. 

4.1. Spectral types 

We take full advantage of the catalog of spectral types 
in the ONC of H97 to determine the temperatures of 
as many members as possible. Matching our catalog 
with the H97 spectroscopy (see Paper I) , we isolated 876 
sources with available spectral type. 

We extend the II97 sample of spectral types using new 
spectroscopic data. The observ ations were ob t ained as 
part of a follow-up to the study ofl Stassun et al.l (|1999f ) of 
the larger Orion Nebula star forming region. The spectra 
were obtained with the WIYN 3.5ni multi-object spec- 
trograph at Kitt Peak National Observatory on UT 2000 
Feb 11. The instrument setup was the same as described 
in H97, providing wavelength coverage over the range 
5000A< A < 9000A with spectral resolution of i? - 1000. 
Exposure times were 2 hours and the typical signal to 
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Figure 9. Same as Figure [S] siiowing tiie difference in modeled 
displacement caused by accretion increasing the accretion temper- 
ature from 7000 K (assumed in the case of Figure [S] solid lines) to 
9000 K (dotted lines). 

noise ratio was ~ 50 per resolution element. Spectral 
types were determined from the spectra using the same 
spectral indices described by H97. In particular, for our 
targets which are principally of K and M type, we found 
the TiO indices at 6760A, 7100A, and 7800A to be most 
useful. We calibrated the indices against observations 
of stars previously classified by H9 7 as w ell as with the 
template library of I Allen fc StromI ()1995D . For 85% of 
the sources a spectral type is assigned with a precision 
of 1 subtype or less. From these observations, we include 
the 65 stars that overlap with the WFI FOV. 

We also add 182 stars classified in Paper I from the 
TiO spectro-photometric index (hereafter [TiO]). These 
sources are limited in the spectral range M1-M6. Earlier 
spectral types are not classified with this method because 
our index correlates with a TiO absorption band present 
only in M-type stars. On the other hand, for late M- 
type stars we were limited by the detection limit of the 
6200A narrow-band filter. Given the small range of Tcff 
spanned, this selection effect in Teg corresponds to a se- 
lection in masses, since the evolutionary tracks are nearly 
vertical in the HRD. 

Thus, spectral classification is available for 1123 out 
of 2621 stars with optical photometry. For 121 of them, 
only /—band ground-based photometry is available, and 
we exclude them from our analysis since the lack of at 
least one color term makes the derivation of the redden- 
ing impossible. The sample of stars with known spectral 
type and at least V and / magnitudes accounts there- 
fore for 1002 objects, 820 of them with slit spectroscopy 
(from both H97 and our new spectroscopic survey) and 
182 M-type stars classified by means of [TiO] index. 

For the stars also present in the H97 catalog, we 
utilize the membership probability collected in that 
work, which is taken from the proper motion survey of 
[Jones fc Walker, (,1988 ) and assigning a 99% probabil- 



ity to the externally photo-evapor ated disks {proplyds^ 
iBallv et al.l 120001 iRicci et al.ll2008[ ) located in the vicin- 
ity of the Trapezium cluster. Among our 1002 stars, 
703 are confirmed members, with a membership proba- 
bility P > 50%, while 54 are non-members (P < 50%). 
The distribution of memberships is bimodal, peaking at 
P = 99% and P = 0%. This implies that our choice to 
define the limit at 50% does not influence significantly 
the selection. In fact selecting only sources with P > 80% 
results in selecting only 12 fewer stars. For 245 the mem- 
bership is unknown, but given the relatively low fraction 
of confirmed non-members, we expect all but 20-30 of this 
stars to be ONC members. This contamination, corre- 
sponding to 2-3% of the sample is too small to bias our 
results; therefore we consider the unknown membership 
stars as bona fide members, taking advantage of the in- 
clusion of 35% more stars in our sample to improve 
the statistical significance of our results. 

4.2. Analysis of completeness 

In paper I we obtained the completeness function of 
our photometric catalog, based on the comparison of the 
luminosity functions (LFs) derived for WFI with that 
of the HST/ACS survey on the ONC, characterized by 
a much fainter detection limit. We refer to this as the 
photometric completeness. Here we compute the com- 
pleteness function of the sample of stars that have at 
least both V and / magnitude in our catalog and and 
available spectral type. We refer to this one as the HRD 
subsample completeness. 

The spectroscopic catalog of H97 is shallower than 
our photometry, with a 50% completeness at Ic — 17.5 
mag, and our extension from new spectroscopy does not 
change significantly this limil^. On the other hand, the 
subsample of stars classified with the [TiO] index is about 
0.5 mag deeper than in H97, to / ^ 18 mag and does not 
cover the bright end of the cluster, populated by earlier 
spectral types (see Paper I). 

The V — I color limit is essentially set by the photo- 
metric completeness at y— band. Our /—band photom- 
etry is much deeper than the y— band: all the sources 
detected in V are also detected in / whereas we detect 
twice as much sources in / than in V. On the other 
hand /?- band is shallower than V , but to provide an es- 
timate of the stellar parameters at least one color (i.e., 
{V — /)) is needed (see Section IT^ . As a consequence, 
the inclusion of a star in the subsample used to derive the 
stellar parameters, is limited by only a) Toff estimates; 
b) y— band photometric detection. 

The 1^— band subsample completeness function is the 
product of the photometric completeness in V and the 
fraction of sources, as a function of V magnitude, that 
have a Tcs estimate. We computed the latter simply 
as a ratio between the LFs in T/— band of the subsam- 
ple and the entire catalog. In this way the subsample 
completeness accounts for both the detection limit of 
the photometry, and the selection effects from the non- 
homogeneity of the sub-sample with Tag estimates. In 
other words, it represents the probability of a star of a 
given y— magnitude to be detected in V and / and to be 

* The new spectra we have presented do not extend the H97 
spectral catalog to fainter magnitudes, and are mostly located out- 
side the FOV of H97, in the peripheral regions of our WFI FOV. 
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Figure 10. The V— band completeness of the sub-sample of stars 
with Tcff estimates (dots). The thick line is a 6th order polynomial 
fit for V > 13.5, while the completeness has been assumed equal 
to a fraction 0.83 for V < 14. 

assigned a spectral type. 

In Figure [TU] we show the HRD subsample complete- 
ness fmiction. We fit a polynomial for V > 13.5, while 
we assume the function to be constant, at 83% complete- 
ness for V < 14: given the high statistical uncertainties 
in the determination of the function. 

We include 50% of the sources at F ~ 20 mag and 
the completeness is zero for V > 22 mag. The com- 
pleteness is lower than 100% for the brightest end of the 
sample, and this is due both to saturated sources (sat- 
urated in / and not present in our catalog) or simply 
because of the lack of spectroscopy. The higher fraction 
atl5<F<18 mag comes from the classification of M- 
type stars from the [TiO] photometric index, capable of 
detecting almost all the sources in this range and there- 
fore providing a more complete spectral typing coverage 
than spectroscopy. 

This completeness is used in the next sections, con- 
verted from being expressed in magnitudes to masses 
and ages, to correct the age distribution and the mass 
function. 

4.3. Spectral type-temperature relations 

In Section HTT] we have collected the spectral types for 
part of the sources present in our photometric catalog. 
These must be converted into Toff, because: 1) for deriv- 
ing intrinsic colors and bolometric corrections (see Sec- 
tion 13. 3|) we rely on synthetic photometry on grids of 
stellar atmosphere models which are computed in terms 
of Teg; 2) also our model for the displacement of the col- 
ors caused by accretion luminosity - which we use in the 
following sections to disentangle this effect for determin- 
ing Ay of the members - is based on the same quantities. 

To derive the effective temperatures from the available 
spectral ty pes we use the s pectra l type vs. temperature 
relation of iLuhman et al.l ()2003fl . This relation, quali- 
tatively intermediate between giants and dwarfs for M- 
type stars, is derived empirically imposing that all the 
four members of the presumably coeval m ultiple system 
GG Tau are located on the same BCAH98 (jBaraffe et al.l 
119981) isochrone (iLuhmanI I1999D . and uses the relation 
from iSchmidt-Kahler I (|1982fr for types < MO . This 
is different from what used in II97, taken largely from 
ICohen fc Kuhii (|1979( ) (and supplemented with appro- 
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Figure 11. Comp arison between the s pectral type-T^ff relations 
used for this work l ILuhman et aL[ )2003l) and the one used on the 
ONC in H97 that was derived in|Cohen & Kuhi (1979). The main 
differences rise with increasing spectral class, and can be several 
hundreds of Kelvin for late-M type stars. 

priat e scales for 0-type s tars (iChlebowski fc Garmanvl 
Il991h and late M-types (|Kirkpatrick et al.l 119931 ) and 
shifted in the K8-M1 range to be smoother). Figure ITT] 
compares the spectral type - temperature relations of 
Luhman with the one of II97. The major differences are 
for the M spectral types, where Luhman predicts higher 
Tcs- 

Using the Luhman relation partially solves a problem 
encountered by H97, the systematic trend of the average 
extinction as a function of temperature with late type 
stars appearing bluer than the color predicted according 
to their spectra type for zero extinction. This effect in 
II97 was partially due to the temperature-intrinsic color 
transformation assumed in that paper, which we already 
improved as described in Section 13.21 but only by us- 
ing the Luhman relation between spectral types and Toff 
is it possible to remove the correlation between mean 
reddening and temperature, characterized by systemat- 
ically lower - even negative - values of Ay moving to- 
wards late-M spectral types. The Luhman scale, predict- 
ing for a given late spectral type a higher temperature 
and therefore bluer intrinsic color, reduces the number 
of sources whose observed colors cannot be matched by 
adding some reddening. A more detailed analysis on how 
the choice of the temperature scale affects the derived ex- 
tinctions is presented in subsection 14.51 

4.4. Derivation of reddening and accretion 

The reddening is defined as the difference between the 
observed color and the intrinsic one, for example the color 
excess E{V - I) = {V - I)ohs - (V - I)o . The intrin- 
sic color {V — I)o is usually derived on the basis of an 
assumed temperature-color relation. In the presence of 
accretion, however, the source colors are also affected by 
Luccr as shown in Figure [51 and since this effect makes 
the colors bluer, partially compensating the reddening 
caused by dust extinction, neglecting the accretion may 
lead to an underestimation of Ay- Furthermore, in our 
analysis of the color-color diagrams we have accounted 
explicitly for surface-gravity effects and for accretion ef- 
fects in all the photometric bands of our photometry. 
This enables us to assess both reddening and accretion 
more realistically than past methods, based on optical 
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Figure 12. Examples of the method applied for the determination of reddening taking into account accretion for three stars in our 
sample. The ONC sample with available T^ff estimates (black dots) is plotted together with our reference model (green line). For each 
plot, the red circle indicates the observed colors of the considered star, and the green squares the predicted photospheric colors (for Ay = 
and Lace = 0) given T^s of the star from its spectral type. The blue curve is the displacement, from the photospheric color, increasing 
Laccr/Ltot from to 1 (see Figure [Sj. The intersection of the reddening direction (dash-dot line) applied to the star, and the accretion 
curve delimits the reddening vector (solid red line). In the three cases we report in the panels the known T^fj and the derived extinction 
and accretion luminosity for the considered star. 



photometry, which assume that only one or the other of 
the effects dominates in the bands in which the observa- 
tions have been obtained. 

We have therefore developed a technique that consis- 
tently takes both Ay and Laccr into account, providing 
at the same time an estimate of the accretion luminosity 
and reddening. Given a star for which T^g is known, our 
reference model gives its original position in the color- 
color diagram for zero reddening and accretion; the mea- 
sured location of the star in the color-color plots will be 
in general displaced because of these effects. But under 
the assumption that they are uncorrelated of each other 
(the differential extinction depends mainly on the posi- 
tion of the star inside the nebula) and since we know how 
to model them (i.e., we know the reddening vector and 
we have a model for the accretion providing the displace- 
ment curve for all possible values of iaccr/^tot), we can 
find the solution that produces the measured displace- 
ment with a unique combination of Ay and iaccr/^tot- 

In practice, for each star in the color-color diagram, 
{B — V) vs. {V — I) we compute the track (like those 
shown in Figure [8]) representing the displacement com- 
puted for increasing Laccr/^tot from to 1 for the exact 
Teff value of the star. Then we found the intersection be- 
tween the reddening vector applied to the position of the 
star and the accretion curve. This provides the exact red- 
dening vector, and its component along the y-axis gives 
the excess E{V — I). The intersection point also gives a 
value of Laccr for the star, according to our model of the 
accretion spectrum. In this way, under the assumption 
that the accretion spectrum we have used is representa- 
tive of all the stars and the reddening law is the same 
for all stars, the values of Laccr and Ay for each star 
are uniquely derived. Some examples of this method are 
plotted in Figure [T2] 



In a few cases this method cannot find a geometrical 
solution; this happens when the star is too blue or too red 
in (B — V), so that the extinction line, applied backwards 
to the star, never crosses the curve computed for its tem- 
perature. Since the photometric errors are higher in B 
than in V or /, we attribute this effect to the inaccuracy 
of the {B — V) color, and therefore we move the point 
along this axis until the closest intersection condition is 
found, respectively the tangency condition of reddening 
line and accretion curve, when the location of the star is 
too on the left. In the diagram, and the starting point 
along the accretion curve where Laccr = when it is too 
on the right. 

For some stars the extinction and accretion determina- 
tion is degenerate. This is the case of Figure [T2b. where 
the reddening vector intersect the accretion curve in 2 
points, either low Ay and low Lace (the solution shown) 
or higher Ay and Lace- In this cases, which involve only 
low-Tefi stars, we always choose the first scenario. 

In a few cases the derived extinction came out nega- 
tive. This happens when the observed {B — V) is higher 
and {V — I) is lower than the locus for zero accretion 
at the star's temperature. For 13% of the sources the B 
magnitude is not available and Ay is therefore derived 
simply from the excess E{V — /). 

In Figure [13] the distribution of Ay is shown, for the 
two choices of Ry. There is a small fraction of stars, as 
mentioned, for which we find a negative extinction. This, 
reasonably accounted for by the uncertainties, is how- 
ever significantly smaller than in HQ?. In these cases we 
consider the star to have zero extinction for the deriva- 
tion of the stellar parameters, as in H97. The median 
of the distribution extinction for the whole sample is at 
Ay ~ 1.5 mag. 

In Figure [14] the distributions of log Laccr/Ltot is 



Optical Photometry of the ONC. II 



13 



150 



P 100 




Figure 13. Distribution of the extinction measured for the ONC 
stars, in terms of Aj in the WFI photometric system, for values of 
Ry = 3.1 and 5.5. The dashed Une shows the median value for the 
extinction for Ry = 3.1, at Ay = 0.92 mag; the dotted line is the 
median for Ry = 5.5 at Ay = 0.97 mag. 
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Figure 14. Distribution of the accretion luminosities derived, for 
the two values of Ry. 



shown, including all the stars for which the value of 
Laccr, derived from our geometrical method, is posi- 
tive and big enough to lead to a displacement in the 
two color diagrams larger than the photometric error. 
We find about 60% of the stars showing an excess at- 
tributed to accretion, with a relative fraction Laccr/.^tot 
peaking at 5-25%. These values are comparable to 
those seen in the more rapidly ac creting Taurus mem- 
bers (jHerczeg fc Hillenbrandl [20081) . Increasing Ry pro- 
duces a decrease in both the fraction of accreting sources 
and the measured Laccr, as a consequence of the steeper 
de- reddening vector in the B,V,I plot. 

We estimate the uncertainties in our determined ex- 
tinction and accretion luminosity using again the ge- 
ometrical method in the 2-color diagrams. Given the 
non-analytical nature of the method, we use a Monte 
Carlo (MC) approach to compute how the uncertainties 
in our data (both photometric errors and uncertainty in 
the spectral types) propagate to the determined values 
of Ay and Laccr- For each star we perturb randomly the 
photometry of each band with their photometric errors, 
assumed to be gaussian. We also consider the uncer- 



tainties in the spectral types from spectroscopy, of ±1 
sub-type, and we assign a temperature within this in- 
terval from a flat distribution. We iterate this approach 
200 times for each star, deriving each time the E{V — 1) 
and log(Laccr/.^tot) values. We then analyze the over- 
all variation of the derived results. Our test shows that, 
on average, <Je{v-i) — 0.12 mag, with the uncertainty 
in the spectral typing being the dominant source of un- 
certainty. Concerning the accretion luminosity fraction 
log(Laccr/.^tot), our MC tcst shows that the uncertainty 
can be up to 1 order of magnitude. We conclude that 
our method is fairly accurate in computing the extinc- 
tion Ay , a most critical parameter for placing the mem- 
bers in the HRD, but provides only a rough estimate of 
the accretion luminosities of the sources. Other methods 
more directly based on accretion indicators (line emis- 
sion, UV-excess,...) should therefore be used for an ac- 
curate analysis of the accretion activity of our sources. 

4.5. Systematic effects in deriving Ay 

In the previous sections we have presented evidences 
that intrinsic colors valid for MS dwarfs are incompatible 
with the observed colors of the ONC, suggesting a gravity 
dependance of the optical color scale for M-type PMS 
stars. This is also predicted by atmosphere modeling, 
though with a limited accuracy. Thus we have calibrated 
the optical colors to best match our ONC data, defining 
our reference model for colors as a function of Tcff . 

We have presented a method to estimate Ay from 3 
photometric bands accounting for possible accretion lu- 
minosity, and we have mentioned that we chose a partic- 
ular spectral type - Toff relation to reduce the number of 
sources showing negative Ay. Here we analyze in more 
detail how all these aspects affect the derived reddening 
distribution. 

We have considered all 8 possible combinations of ei- 
ther considering or not considering the 3 described im- 
provements, which i s: a) converting the spectral types 
into Tpff using the iLuhman et aLl (|2003f) scale or the 
ICohen fc Kilhil (11979^ one, used, e.g., in H97; b) using 
the intrinsic colors in B, V, I band from our reference 
model or from the log g values of the ZAMS; c) deriving 
Ay using the from BVI disentangling accretion as de- 
scribed in Section [4.41 or simply neglecting the accretion 
and the i3-band information and basing solely on the ex- 
cess E{V — I). In the case of the ZAMS model, we have 
first computed the displacements in the 2-color diagrams 
caused by accretion luminosity, in a similar way as for 
our reference model, shown in Figure [5] 

The results are presented in figure [13 as a function 
if Toff to highlight systematic trends. For guidance, the 
top-left panel shows the results we use, and the bottom- 
right one reproduces qualitatively the findings using the 
same methods and relations as in II97, with a systemat- 
ical lower Ay towards the latest M-type stars. 

It is evident that using the Luhman temperature scale 
for M-type stars instead of the Cohen & Kuhi relation 
leads to systematically higher Ay, compensating very 
well the problem of deriving negative reddenings. 

The differences we observed in the Ay values changing 
the color scales from ZAMS to our empirical reference 
model are more modest, and Figure [15] alone does not 
allow to clearly justify the choice of one or another color 
scale. We remind the reader that the intrinsic colors from 
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Figure 15. Ay vs. log(Tcfl) plots for the ONC members, for all the combinations of changing: a) Tcffvs. spectral type scale - left panels 
using the Luhman scale, right panels using the Cohen & Kuhi scale; b) intrinsic colors: - rows 1 and 2 using our empirically calibrated 
reference model, rows 3 and 4 using the result of synthetic photometry constraining logg from the Siess ZAMS; c) method to derive Ay'- 
- row 1 and 3 using our geometrical method to disentangle Ay from Laccr, row 2 and 4 neglecting B magnitude and deriving Ay from 
E{V-I). 
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the ZAMS model do not reproduce the observed color- 
color diagrams, (see Figure [2]), and therefore can not be 
considered representative of the observed data. In partic- 
ular, assuming the dwarfs colors leads to systematically 
higher accretion luminosities for cool members, predict- 
ing almost all M-type stars to have a positive value of 

tot- 

Finally, if we do not consider the color offsets due to on- 
going accretion (Figure [151 second and fourth row), the 
predicted extinctions are underestimated. This is clear 
from figure [5] a positive value of Laccr leads to lower 
intrinsic colors (both B — I and V — I), increasing the 
measured Ay- 

5. THE H-R DIAGRAM OF THE ONC 

The Hertzsprung-Russell diagram, representing the 
total luminosity emanating from the stellar photo- 
sphere versus the photospheric effective temperature, 
is the physical counterpart to the observational color- 
magnitude diagram. The transformation between the 
two requires the knowledge of a number of quantities. In 
our case, having at our disposal Toff, we convert the ob- 
served magnitudes into total luminosities correcting for 
accretion and extinction (Section 14. 4p . using our com- 
puted bolometric corrections for the WFI photometric 
system (Section 13. 3p . and considering the solar colors 
and bolometric magnitude. We therefore compute the 
luminosity as follows: 

log (^-^^ = 0.4 ■ [Mboi,e -Mboi,.]{7) 

= 0.4 ■ [M/^pj,0 — M/„pj,^, -fBC/„pj,0 — i3C/„pi,^] 

= 0.4 • [Mv,Q -{V- /wFi)o - %Fi + AIacc + 

Ai„^, + BCi„^,iT) + DM] 

where /wfi and Vyvpi are magnitudes referred to the 
WFI system, while V is Vjohnson- The term A/accr (nega- 
tive) is the / band excess due to the accretion luminosity 
derived as shown in ij4.41 directly computed from syn- 
thetic photometry on the accretion spectrum normalized 
to the measured value of iacc of each star, for the 2/3 
of the stars showing accretion in their colors. A/accr has 
a mean value of— 0.16±0.1 mag, and is clearly zero for 
the rest of the sample. According to our convention for 
the bolometric correction, we have BC/„pj,Q= 0. For 
the solar values we used Mvq = Afy(johnson)0 — 4.83 
(jBinnev fc Merrifieldl Il998f ) and by means of synthetic 
photometry on our solar spectrum we find (Vjohnson — 
-'^WFi)© — 0.80. Regarding the distance of the ONC, re- 
cent investigations based on accurate trigonometr ic par- 
allax measu rements (Hirota ct al. 2007; Sandstro m et al.l 
I2007t iMenten et al.ll2007l) lead to a revision of the value 
d ~ 470 pc adopted by H97. Amongst t hese, we used 
the m ost precise estimate d = 414 ± 7 pc (|Menten et al.l 
[2001 . based on a recent VLBA trigonometric parallax, 
corresponding to a distance modulus DM= 8.085. 

As mentioned in Section [4l we have excluded from our 
sample all the confirmed non-members, and we expect 
that the remaining contamination from foreground and 
background sources with unknown membership is rela- 
tively low (^2 - 3% of our sample). 

The resulting HRDs, for the standard Ry = 3.1 red- 
dening law, are shown in Figure 1161 We mention that 



assuming the higher value of Ry ~ 5.5 leads to pre- 
dicted higher luminosities because, for a given color ex- 
cess E{V — I) the extinction Aj is higher (see Table [2]). 
The upper and lower panels in Figure [TBI correspond to 
different PMS evolutionary models; top) the models by 
Siess, with ages from 10^ years to the ZAMS; bottom) 
the models of lPalla fc Stahleii (|1999D (ages from 10^ to 
10* years), which include the deuterium birthline. Our 
sample reaches masses as low as O.IM©, extending to the 
intermediate-mass range up to several solar masses. 

Figure [161 highlights the differences between the two 
families of evolutionary models in the comparison with 
our data. Whereas both are generally parallel to the 
median of the observed sequence, Siess predicts that 
isochrones of a given age are located at higher luminosi- 
ties, which implies t hat the ages predicted for the cluster 
members are older. iPalla fc Stahleii (|1999[ ). on the other 
hand, predict lower luminosities. 

There are a few anomalous stars located in regions of 
the HRD not covered by the PMS isochrones. Stars lo- 
cated below the ZAMS can be background objects with 
no membership information. However, there are also 
some u nderluminous confirmed mem bers. As explained, 
e.g., in iKraus fc Hillenbrand! (|2009f ). this effect seen in 
other PMS systems is due to the presence of dense cir- 
cumstellar material (e.g., an edge-on circumstellar disks) 
which may both a) occult partially the photosphere; and 
b) scattering a significant amount of stellar flux - more 
efficiently for short optical wavelengths - along the line 
of sight. This hypothesis has rec eived recent support the 
study of iGuarcello et all (|2010l ). All the sources below 
the Siess ZAMS or the 100 Myr Palla isochrones are ex- 
cluded from our further analysis. 

We assess the overall uncertainty, from the observa- 
tional errors, in the position of the stars in the HRD, 
using a Monte Carlo approach similar to that adopted in 
Section 14.41 We consider the uncertainty in the spectral 
type and the photometric errors, apply these effects ran- 
domly a number of times to each star, deriving the final 
position in the HRD as described above. For each source 
the methods produces a "cloud" of points in the H-R 
diagram which represent the real scatter associated with 
the star. Figure [TTl shows the result of an iteration. Each 
star is associated to a stripe of points, whose inclination 
and shape in general depend on the location of the star in 
the 2-color diagrams and its temperature. To understand 
better the general behavior of the resulting uncertainties, 
we divide stars in bins of log(Tcff) and compute the av- 
erage spread in both T^e and logL with respect to the 
actual position of the stars of Figure [161 This is repre- 
sented with black points and error bars in Figure 1171 It 
is evident that the uncertainty in T^s, which comes di- 
rectly from the uncertainty in assigning spectral types, 
is usually much higher than the spread in log(L), even if 
this is infiuenced by the former (as mentioned in Section 
14.41 the errors in Tcff usually contribute more than the 
photometric errors in the derivation of Aj, which affects 
log(i). 

Concerning the luminosities, however, we must stress 
that there are several sources of uncertainties beyond the 
measurement errors, which involve both the data anal- 
ysis and the physical nature of the members. In fact, 
our determination of Ay and Laccr from the BVI pho- 
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Figure 16. Hertzsprung-Russel diagram for the ONC. Isochrones and tracks from lSTess et al] II200CII ') (upper panel) and lPalla fc Stahleil 
(|T999) (lower panel) are overplotted, for masses from O.IMq to 3.5Mq and ages from 10^ years to the ZAMS in the Siess case, and from the 
birthline to 10* Myr for the lPalla fc Stahleil I I1999I ) models. Blue stars are spectral types from H97, excluding sources with a membership 
probability < 50%. Red circles indicate the additional sources we classify from our spectroscopic survey. Green crosses show the location 
of the M-type source, for which the temperature is determined from the [TiO] spectrophotometric index (see Paper I). 



tometry relies on the assumption of a unique accretion 
spectrum in the modehng of the veihng. The SED of the 
accretion spectrum - and therefore its broad band colors 
- may actually be different from star to star. As men- 
tioned in Section 13.51 these variations should not affect 
much the extinction, and so Ltot- We do not have at 
our disposal a representative sample of calibrated accre- 
tion spectra for the ONC, mostly due the bright nebular 
background against which a flux-calibrated excess would 
have to be measured relative to the stellar continuum. 
Therefore we are unable to characterize explicitly the 
error distribution on the luminosity caused by the ac- 
cretion spectrum modeling. Other sources of apparent 
luminosity spread are, moreover, circumstellar emission, 
scattered light, unresolved companions, variability and 
accretion history. As a consequence, the errors on Ltot 
shown in Figure [T7] do not represent the maximum uncer- 
tainty, but the minimum error based on the uncertainties 
of our data that are known. Only a quantitative knowl- 



edge of all the sources of luminosity spread could disen- 
tangle, therefore, apparent luminosity spread from true 
age spread. This problem, critical for the understanding 
of age distribution in the star formati on histor y of young 
stellar systems (e.g. , Hillenbr and et al . 2008; H illenbrandl 
I2009t IBaraffe et all 12009.) will not be discussed in detail 
in this work. 

We can compare our derived luminosities with the val- 
ues found by H97, for the sources present in both cat- 
alogs. The results, as a function of Toff, are shown in 
Figure [JHl Besides a modest scattering in the result, 
most probably due to stellar variability, a definite trend 
is evident: for log Toff > 3.6, corresponding roughly to 
K and earlier spectral types, we estimate a lower bolo- 
metric luminosity. This is in part caused by the lower 
distance we adopted for the ONC; however this accounts 
for only 0.11 dex of logL. The changes we applied to 
the intrinsic colors and to the temperature scale do not 
affect critically this range of spectral types; therefore. 
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Figure 17. Representation of the method we used for determining 
the average uncertainties that affects our determination of the H-Il 
diagram. For each star we propagate the uncertainties simulating a 
"cloud" of determined positions in the diagram. The black points 
with error bars represent the average determined uncertainty in 
logL and logT in 12 bins of surface temperatures. 
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Figure 18. Difference between bolometric luminosities derived 
from our calculations and from Hillenbrand (1997), as a function 
of Teff. The systematic trend is evident: for early type stars we 
derive a lower value of L, while for M-type stars our luminosities 
trend towards increasingly higher values towards lower tempera- 
tures, relative to H97. A number of differences in the analyses 
contribute to this trend. 

this systematic effect is largely due to our bolometric 
corrections from /—band magnitudes. A different type 
of discrepancy appears for cold stars, where our analysis 
leads to higher luminosity than in H97. This is mostly 
due to the higher extinctions compared to H97 that we 
measure in this range. Since the Luhman temperature 
scale we adopt predicts higher Toff in this interval, our 
HRD is shifted both to higher log L and higher log T for 
M-type stars compared to II97, and continue to follow 
fairly well the isochronal slopes of evolutionary models. 

5.1. Completeness of the HRD 

In section 14.21 the completeness of our spectro- 
photometric catalog is derived as a function of V mag, 
i.e. for the CMD. This function, however, can not be 
directly translated in the HRD, because of differential 
extinction and non-linearity of intrinsic colors and bolo- 
metric corrections. Therefore we developed a statistical 
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Figure 19. Simulated completeness in the HRD 

simulation that, for each point of the CMD, determines 
the probability for such star to be observed and present 
in our spectral type sample. 

We proceed as follows. We consider a uniform, dense, 
cartesian grid of points in the HRD. For each of these 
points we find its counterpart in the CMD for Ay = 0, 
converting Tcff and logL into {V — I)o and Iq using the 
intrinsic colors of our reference model and applying back- 
wards the bolometric correction and the distance modu- 
lus. In reality, a star having the considered Toff and log L 
of the grid point is not necessarily observed in this posi- 
tion of the CMD, but reddened along the extinction di- 
rection of an amount that follows the distribution of Ay 
inside the nebula. We consider as a representative distri- 
bution of Ay the distribution we determined in section 
[44] limited for the brightest half (in logL) of the popula- 
tion. This is because the reddening for intrinsically faint 
stars in our sample is biased towards low values, since 
such sources with high Ay fall below our detection limit. 
For each point in the transformed grid in the CMD we 
utilize a Monte Carlo approach for applying the extinc- 
tion, creating 1000 simulated reddened stars, with values 
of Ay drawn from this distribution. Therefore, this set 
of aligned, simulated stars correspond to the all the po- 
sitions in the CMD where a star with the considered Toff 
and log L can be located, following the true extinction 
distribution. We assign the completeness computed in 
Section |4?2] to all the 1000 simulated stars according to 
their position in the CMD. The average of the 1000 values 
of completeness correspond statistically to the real prob- 
ability, for a source with the considered Tcff and log L, to 
be detected in our catalog and assigned a spectral type. 
Iterating the same method for all the points in the HRD, 
we populate a 2D distribution of completeness in this 
plane. 

The result is shown in Figure 1191 Our results show 
that the completeness decreases towards both low L/Lq 
and low Toff. 

5.2. Distribution in the mass-age plane 
Using the n iodels of iSiess et all (j2000( ) and 



IPalla fc Stah"ieJ (fT999h . we assign an age and mass 
to each star. To this purpose, we compute 2D interpo- 
lated surfaces with age and mass given in function of 
log(Teff),log(L/T0). In Table [3] we present our results 
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for the stars of our sample. 

We analyze our masses and ages, and study their dis- 
tribution in the mass-age diagram; this allows us also to 
investigate possible correlations between the two quanti- 
ties. 

For this purpose, we produce a number density map 
in a \og{age) vs. log(Af) plot, counting the stars present 
into uniform bins of width 0.25 dex in log{age) and 0.1 
dex in log(M). A boxcar smoothing with a kernel size of 
1.3 times the original grid size is applied to smooth out 
local variations due to statistical noise. 

The result is shown in Figure 1201 upper panels, with 
density contours highlighted. The outermost solid con- 
tour corresponds to 1 star per bin, the second 5, and the 
third 9, corresponding to poisson uncertainty of 100%, 
45%, and 33%. Evidence of correlations should be in- 
vestigated in the inner part of the map, where the ma- 
jority of the objects are positioned. Here it is evident 
that whereas Siess models do not show a global trend, in 
the Palla & Stabler case the low-mass stars turn out to 
be systematically younger than the intermediate masses, 
with an average age changing by more that one dex along 
our mass range. H97 reported a similar inconsistency 
in their analysis, based on t he evolutionary models of 
iD'Antona fc Mazzitelll (|1994[ ). The fact that the age- 
mass correlation is found using Palla & Stabler mod- 
els while not seen in our results from Siess isochrones 
could imply that the latter should be considered in bet- 
ter agreement with our dat a. 

However, as discussed in IPalla &: Stahleil (jl999D . the 
completeness of the sample can potentially bias the mass- 
age relation. This is especially true in the very-low mass 
regime, as seen in Figure [191 as older PMS stars have 
lower luminosity than younger objects and therefore are 
more likely to be excluded from a stellar sample due to 
sensitivity limits. We therefore investigate if the incom- 
pleteness of our stellar sample could be responsible for 
the observed mass-age relation in order to correct the 
derived results. 

We compute the completeness in the mass-age plane 
following an analogous method as for the HRD (Section 
15. 1|) , but considering instead of a uniform grid of logTcg 
and logL values, a cartesian grid of masses and ages 
in the HRD, identical to the one used for Figure [^Db.b. 
The result is now model-dependent, since a given mass- 
age point is located in different positions of the HRD 
according to the two sets of evolutionary models. 

The result is shown in Figure [20l middle panels, inter- 
polated and smoothed in the same way as for the density 
maps. Results confirm that completeness decreases to- 
ward lower masses and higher ages, in a qualitatively sim- 
ilar fashion to the observed mass-age correlation found 
from our data with Palla & Stabler isochrones. Differ- 
ences between the results obtained with the two sets of 
evolutionary models are due to the different shape of 
tracks and isochrones, especially in the low-mass range 
(see the HRDs of Figure [T6l) . The cutoff for the youngest 
ages in the completeness derived using Palla & Stabler 
models is due to the lack of the early PMS phases in this 
family of models. 

The "normalized" mass- age map, computed by divid- 
ing our preliminary results of Figure [^Ob.b by the com- 
pleteness map, is shown in the bottom panels of the same 
Figure. While the shape of the distributions in the mass- 
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Figure 20. Upper panels: density map in the plane log(age) vs. log(mass) according to our results using ISiess et al.l 1 I2OOOD and 
IPalla &: Stahlej I I1999I ) evolutionary models. Stars are counted into a grid of uniform bins spaced 0.1 in log(M) and 0.25 in \o^(age), and 
the derived map has been oversampled, interpolated and smoothed with a kernel sized 0.13 X 0.32. The color scale reported is in units of 
number of stars per unit of log(M) and per unit of log{age). It is evident that in the Siess case the age distribution is uniform (besides a 
broad age spread at all the masses), but for the Palla &: Stahler models there is an evident correlation between age and mass, predicting 
younger ages for very low mass stars and older for intermediate masses. Middle panels: contour plots of the completeness function in 
the mass-age plane, computed for the evolutionary models of ISiess et al.l 120001 ) and IPalla &: Stabler! 119991) . and for Ry =3.1, from our 
simulation. The apparent vanishing completeness in the bottom-left corner of the Palla &; Stahler case is due to the lack of information of 
these evolutionary tracks for young ages and very low masses, but this does not affect our results given that our sample includes no stars 
in this range. Bottom panels: the density maps of the upper panels normalized dividing it by the completeness. 
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age plane and the average age for the ONC do not change 
by a large quantity, due to the smooth variation of com- 
pleteness against mass and age, we find that the cor- 
relation between mass and age for the Palla & Stabler 
models reduces significantly, and below IMq the trend 
of age with mass is less evident, while, on the contrary, 
for the Siess models the completeness correction suggests 
an anti-correlation between ages and masses, hinting for 
higher predicted ages for low-mass stars. 

In the high-mass end of these plots (M > IM0) there 
is still a tendency toward a higher measured average age, 
but we do not consider this as a solid evidence for a mea- 
sured mass-age correlation. This for two main reasons: 
on the one hand, intermediate stars of young ages are not 
covered by the evolutionary models of Palla & Stabler, 
and on the other hand for low-mass and high-age fraction 
of the population the completeness we determined is too 
low to be correctable using our method. In conclusion, 
we do not find evidences of an age trend with respect to 
stellar mass in the ONC, both using Siess and Palla & 
Stabler evolutionary models. 

The apparent over-density evident in the completeness- 
corrected map at logM — O.SM© for log(age) > 7 is not 
relevant, being solely due to statistical noise. 

5.3. The age distribution 

In Figure [21] we present the measured age distribu- 
tion and the same after correcting for the completeness 
C(t), for the two families of PMS models, as we now 
describe. From the 2D completeness functions in the 
mass-age plane, we can now derive the overall complete- 
ness in ages. If C(M, r) is the completeness function of 
mass and age, shown in Figure [20b.f. the completeness in 
ages C(r) is just its projection normalized on the density 
in the mass-age plane: 



C{t) 



J C{M, t) ■ n{M, t) dM 
/ n(A/, r) dM 

^ / n,amplc(M,T) AM 

J n{M, t) dM 



(8) 



where n(M, r) is the completeness-corrected density in 
the mass- age plane fFigure BDb.f ) and nsampio(-M, r) the 
measured one (Figure [20k. b) . 

The completeness correction does not infiuence signif- 
icantly the average ages and the measured age spread; 
but increases t h e old age tail of the distributions. In 
the iSiess et all (|2000l ) case the age distribution peaks 
at ^ 3 Myr, w i th a broadening of 0.3 dex while with 
[Palla fc Stabled ([1999[ ) the peak is at about 2 Myr with 
a higher broadening of 0.4 dex. However, the broader 
distribution of the latter in logarithmic scale vanishes in 
linear units, because of the younger average age. Within 
Icr the age of the population spans from 2.5 to 5 Myr for 
Siess, and from 1.22 to 3.22 Myr for Paha & Stabler. 
Increasing the slope of the reddening law produces a 
modest shift towards younger ages (< 0.1 dex)because 
of the higher luminosities derived. As mentioned in the 
previous section, the observed age spread does not co- 
incide with a real star formation history, because of the 
unknown uncertainties in the stellar luminosities. There- 
fore the distributions shown in Figure [21] represents the 
combined effect of the real age spread and the apparent 
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Figure 21. Age distribution for the ONC population. Top panel: 
ages derived using the Siess et al. (2000) evolutionary models; bot- 
tom panel: the same using Palla & St ahlerl (ll999i V The filled his- 
tograms show the distributions of ages in our sample. The open 
histograms are the same after correcting for incompleteness, and 
the statistical uncertainties are shown for each age bin. 



spread originated by other effects that scatter the mea- 
sured luminosities. The latter is mainly due to stellar 
variability and scattered light from circumstellar mate- 
rial, and is probably responsible for the 7-9% (respec- 
tively for Palla & Stabler and Siess models) of stars with 
a measured age exceeding 10 Myr, percentage that in- 
creases of ^ 50% when correcting for completeness. This 
fraction of old age sources is about 3-5 times higher than 
the upper limit for the residual contamination from non- 
members (see Section [3|). 

We stress that the sources of apparent luminosity scat- 
tering we have mentioned may be responsible for the tails 
of the age distribution, but we do not expect the overall 
width of the age distribution to be significantly overes- 
timated. In particular, our a ge spread agree s with the 
one measured in the ONC bv [ Jeffrie^ (|2007| ) by means 
of a statistical analysis of geometrically estimated radii - 
therefore immune of the effects of hidden binarity, vari- 
ability, extinction uncertainties. This suggests that a real 
age spread is the main cause of the observed luminosity 
spread in the H-RD and that we are not significantly 
overestimated the width of the age distribution. 
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Figure 22. Completeness as a function of stellar mass computed 
for both evolutionary models. 

5.4. The mass function 

The mass function (MF) of a system is defined as the 
number distribution of stars as a function of mass. While 
a general observational approach measures the so-called 
present day mass function (PDMF), a cornerstone for 
understanding how stars form is the Initial Mass Func- 
tion (IMF), which is the mass distribution according to 
which stars are born. For a pre-main sequence popula- 
tion younger than the time required to produce a sig- 
nificant mass segregation, or to evolve the most massive 
stars into post-MS phase, the observed mass function co- 
incides with the IMF. 

The derivation of the IMF from observations always 
presents several sources of uncertainty, the main ones 
typically being unresolved companions and completeness. 
We can asses the completeness of our sample as a func- 
tion of mass from the 2D completeness in an equivalent 
way as for the age completeness (eq. [8]), projecting it on 
the mass scale: 



C(M) 



J »samplc(M,T) dr 

/ n{M, r)dr 



(9) 



The results, computed for both evolutionary models, 
are presented in Figure [221 In the intermediate-mass 
regime both evolutionary models provide similar com- 
pleteness, whereas for lower masses the Palla & Stabler 
models lead to a systematically lower completeness. 

To derive the IMF of the ONC we count the number 
of stars in equally spaced bins 0.1 dex wide in log M and 
normalize the result to the bin width in order to express 
the IMF in terms of ^(Af) instead of ^(logM). We cor- 
rect the star counts for completeness, interpolating the 
function shown in Figure [52] to the mass values of the 
IMF. 

The two panels of Figure [23] show the IMF derived 
for both Siess and Palla & Stabler evolutionary mod- 
els. In the same figure, the H97 IMF is overplotted, 
showing that it follows our distribution f or intermediate - 
and low-mass stars (M > 0.4Mq). The IKroupal (|2001[ ) 
and ISalpeteri (jlOSSl) IMF are shown as well, evidenc- 
ing the differences with our mass distributions. Specif- 
ically, Kroupa under-produces the observed low-mass 
population and does not exhibit the observed turn-over. 
Salpeter over-produces the observed number of low mass 
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stars with the effect perhaps even larger than indicated 
due to the normalization in the mid-masses rather than 
the high mass end. Table [4| provides the value for £,{M) 
and the associated uncertainty for the two cases. 

As we mentioned in Section [5.3[ the luminosity spread 
we measure can be slightly overestimated due to phys- 
ical and observational biases such as stellar variability, 
hidden binarity, scattered light, and uncertainties in the 
positioning of the members in the H-RD. However, since 
the evolutionary tracks in the H-RD are nearly vertical 
in the mass range we investigate, and given that, as men- 
tioned, a true age spread is probably the dominant source 
of luminosity spread, we do not expect that our IMF is 
significantly affected by these biases. 

The two sets of models produce significant differences 
in the IMF for low-mass stars, log(Af) < O.SMq. For 
both models the IMF shows a fiattening at low masses. 
This feature is real, given our accurate assessment of 
completeness. However, while for Palla & Stahler the 
change of slope from the general power-law is relatively 
modest and in agreement with the Kroupa mass distri- 
bution, the Siess model leads to a clear turn-over be- 
low 0.2Mq, compensated by an apparent overabundance 
of stars in the range O.2M0 < M < 0.3Mq. Similar 
differences in the shape of the IMF below 0.2Mq, due 
to differences between theoretical models, temperature 
scales and bolomctric corrections, were also reported by 
Hillenbrand & Carpenter (2000). 

We tested also the changes when Ry is increased (not 
shown in the Figure), finding that this does not intro- 
duce significant modifications in the shape of the IMF. 
This is not surprising, considering that the reddening 
law affects the luminosity of stars, shifting the point in 
the H-R diagram along the y— axis in proportion to the 
color excesses E{V — I), whereas for most of the stars in 
this temperature and luminosity range the evolutionary 
tracks are nearly vertical. 

Wit h respect to the IMF presented by iMuench et al.l 
(j2002f ). derived from a Monte Carlo best- fit to the ob- 
served fiT— band luminosity function, we confirm the 
power-la w increas e with slope similar to the x — —1.35 
value of ISalpeted (flOSl . for masses larger than M ~ 
0.6 Mq. However, the break they found at M ~ 0.6 M© 
occurs in our case at lower masses, M ~ 0.3 Mq or less. 
This is consistent with the fact that there ONC members 
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Figure 23. Complet eness-correc ted mass functions derived for the ONC (filled squares) according to the models of lSiess et al.l l|2000l ) (left 
panel) and , Palla fc Stahleil I I1999I ) (right panel). The open circles are the same without the completeness correction. Units are number of 
stars per solar mass (i.e., the counts have been normalized to the v^idth of each bin in units sola r mass es). The mass function for the ONC 
d erived in H97 is superimposed to highlight the differences in the low mass range. The lSalpeteil l|1955l ) single power law IMF and the IMF 
of lKroupal I I2001I ) are shown as well. 



are known to present band excess due to circumstel- 
lar material, which biases the fluxes towards brighter val- 
ues which are erroneously interpreted as higher masses 
than the underlying stars truly possess. 

Thanks to the highest number of sources present in 
the subsample used to derive the IMF, as well as our 
completeness correction, we are able to constrain the 
position of the IMF peak at low stellar masses with a 
higher confide nce than previous works. As discussed in 
iBricenol (pOOl . the presence of a peak in the stellar dis- 
tribution at spectral types ~M3 (corresponding roughly 
with masses M r^ O.SM©) is measured for all the young 
populations of the Orion complex, suggesting a common 
origin for all these populations. On the contrary, other 
star forming regions such as Taurus and IC348 present 
IMFs increasing down to the hydrogen burning limit, and 
in the case of Taurus, showing evidences of stellar over- 
abundance at M r^ 0.6-0.8Mq. According to our results, 
we confirm the presence of a peak in the very-low mass 
regime, but this is located lower masses, closer to 0.2Mq. 

In conclusion, our analysis shows that especially at the 
lowest stellar masses, and presumably also in the entire 
brown dwarf regime, the actual shape of the IMF de- 
pends strongly on the model assumptions. Our current 
understanding of the pre-main-sequence stellar evolution 
suffers from the lack of a general consensus on the best 
model to adopt. By combining accurate and simultane- 
ous multicolor photometry with spectroscopy it is possi- 
ble to derive critical information to constrain the models, 
as we have tried to show in this paper. An extension of 
this analysis in the substellar mass range will be enabled 
by the study of the HST ACS photometry {Robberto et 
al., in preparation) obtained for the HST Treasury Pro- 
gram on the Orion Nebula Cluster, to which this ground- 
based observational effort also belongs. 

6. SUMMARY 



In this work we perform an analysis of the simulta- 
neous multi-band op tical p hotometry from WFI obser- 
vations presented in iDa Rio et al.l (|2009b( ). deriving a 
new H-R diagram of the ONC population in the low- 
and intermediate-mass range and the corresponding mass 
function. 

From the comparison of the observed color-color dia- 
grams obtained from B,V,I and 6200A narrow-band fil- 
ters, with results of synthetic photometry, we find evi- 
dences that the intrinsic colors for the ONC differ from 
those of main-sequence dwarfs in the M-type regime, but 
they are shifted towards bluer colors, in a similar way as 
predicted by atmosphere models assuming logg values 
from a PMS isochrone. Nevertheless, the atmosphere 
models we use are still unable to match accurately the 
photospheric emission of the ONC; we therefore calibrate 
empirically the intrinsic color scale based on our ONC 
photometry, limiting to the known members with modest 
extinction and no accretion detected from Ha photome- 
try 

The color-color diagrams of the population also sug- 
gests that the typical Galactic value of the reddening 
parameter Ry = 3.1 is more compatible with our data 
of the ONC than a higher value of Ry = 5.5. 

Our derivation of intrinsic colors, in particular, allows 
us to improve the accuracy of the derived Ay for the 
ONC members, a critical step for a precise positioning 
of the stars in the HRD. We find evidence of i3-band ex- 
cess for a fraction of sources in the color-color diagrams, 
which we attribute to stellar accretion. We develop a 
method that, using 3 bands simultaneously, allows us to 
disentangle the effects on the broad-band colors of red- 
dening and accretion luminosity. We test the accuracy 
of this method finding that while it leads to a robust and 
more unbiased determination of Ay , it is still insufficient 
to provide a strong constrain on the accretion luminosity. 

The derivation of all the stellar parameters still needs 
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an ind ependent estimate o f Teg . We use the spectral cat- 
alog of lHillenbrandl (|1997[ ). as well as new spectral types 
that we present here, and a sample M-type stars classi- 
fied by means of narrow-band photometry (Paper I) . This 
allows us to have at disposal the most complete sample 
of stars with known spectral type in the ONC. We study 
how the choice of the spectral type-Teff relation is critical 
to remove syste matic trends in Ay as a function of Tag. 
We find the the lLuhman et a l. (2000) temperature scale 
improves over the lCohen fc~Kuhii (|197l scale. 

We derive the bolometric luminosities for ~ 1000 mem- 
bers and position them the HRD, down to the hydrogen- 
burning limit. The revised lower distance we used of 
414 pc compared to previous works leads to systemat- 
ically lower luminosities; this is compensated for in M- 
type stars, where we predict higher L, caused by higher 
measured extinctions for late-type stars. This is due to 
our improvements in the intrinsic colors and our method 
to account also for accretion effects in the observed col- 
ors. 

We compute the 2D completeness function of our sam- 
ple of members in the HRD, including both photometric 
detection and the spectral subsamp le. Using the evolu- 
tionar y models of iSiess et al.l (|2000[ ) and lPalla &: Stabler! 
we assign ages and masses to all the members. 

We extract the age distribution from the different 
choice of models. In the case of Siess the average ages 
turn out to be higher (~ 3 Myr) than Palla & Stabler 
(~ 2 Myr) with a slightly lower spread (0.3 dex vs 
0.4 dex). We find a mass-age correlation if one uses 
the Palla & Stabler models, with high-mass stars ap- 
pearing older than low-mass stars. This trend, however, 
disappears when a careful statistical analysis of the com- 
pleteness of our stellar sample in the mass-age plane is 
performed. 

We derive the IMF of our sample, corrected for detec- 
tion incompleteness down to the hydrogen burning limit; 
the IMF turns out to be model-dependant in the very- 
low mass range, with Siess predicting fewer low-mass 
stars and a turn-over at M ~ 0.2Mo while using Palla 
& Stabler a more modest flattening is present. 
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